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In this paper, we study the role of collective vibrational motion in the phenomenon of electronic energy 
transfer (EET) along a chain of coupled electronic dipoles with varying excitation frequencies. Previous exper- 
imental work on EET in conjugated polymer samples has suggested that the common structural framework of 
the macromolecule introduces correlations in the energy gap fluctuations which cause coherent EET. Inspired 
by these results, we present a simple model in which a driven nanomechanical resonator mode modulates the 
excitation energy of coupled quantum dots and find that this can indeed lead to an enhancement in the trans- 
port of excitations across the quantum network. Disorder of the on-site energies is a key requirement for this 
to occur. We also show that in this solid state system phase information is partially retained in the transfer 
process, as experimentally demonstrated in conjugated polymer samples. Consequently, this mechanism of vi- 
bration enhanced quantum transport might find applications in quantum information transfer of qubit states or 
entanglement. 

PACS numbers: 05.60.Gg,78.67.Hc,07.10.Cm,71.35.-y 



INTRODUCTION 

The rapid and efficient transfer of optical excitation energy 
through a network of coupled dipoles has recently become the 
subject of intense study in a number of natural and artificial 
systems. A typical naturally occurring example is the phe- 
nomenon of electronic energy transfer (EET). Many natural 
photosynthetic proteins and artificial polymers possess light- 
absorbing molecular dipoles (chromophores) which interact 
with each other causing the transfer of energy from an ini- 
tially excited donor to another acceptor chromophore . Ar- 
tificial systems of coupled semiconductor quantum dots can 
exhibit similar energy transfer|0,[3tl- The proper understand- 
ing of EET is important for the study of complex quantum 
systems and may lead to the design of more efficient devices 
for energy storage and conversion. 

A number of papers in which the authors have tried to un- 
derstand the basic mechanisms involved in efficient EET have 
lately appeared in the literature The role of quantum 

coherence in the efficiency of excitation transfer in coupled 
systems is studied in 101 ■ The role of decoherence processes 
on energy transfer within photosynthetic systems were inves- 
tigated in |5, 6j where noise-assisted quantum transport is in- 
troduced. The idea that transport of excitations across chro- 
mophoric networks can be enhanced by local dephasing noise 
was further explored in ITjrtSO where it is applied to the de- 
scription of excitation transport across the Fenna-Matthew- 
Olson (FMO) complex 11311 . This system acts like a 'quantum 
wire' connecting excitations in an antenna complex, where 
light is absorbed and an exciton created, to another com- 
plex that contains the reaction center where the exciton en- 



ergy is used to initiate photo-synthetic chemical processes. 
These studies suggest how nature might use noise-assisted 
processes to increase the efficiency in the performance of cer- 
tain tasks. These concepts have recently been extended to 
the non-Markovian regime|fT5]. An important experiment on 
EET has recently reported an investigation of the intermedi- 
ate coupling regime lITsll . This regime interpolates between 
the classical incoherent hopping (Forster) mechanism and a 
coherent (wavelike) quantum mechanism. The intermediate 
coupling regime is particularly interesting because both popu- 
lation and phase information is supposed to be partially trans- 
ferred through space. Needless to say, such a type of coherent 
EET would be very interesting for applications in the area of 
quantum information. 



In IllSn . coherent intrachain EET in conjugated polymer 
samples are studied using specially designed experimental 
techniques (two-time anisotropy decay (TTAD) and 2D pho- 
ton echo experiments) that enable the measurement of coher- 
ent EET directly and provide information about coherences. 
Their results show evidence of coherent EET in the interme- 
diate regime (phase information is partially preserved). Pro- 
vided that fluctuations of the electronic transition frequencies 
of the donor and acceptor chromophores are correlated, coher- 
ent EET can be observed in two-time anisotropy decay mea- 
surements ifTsll . On the contrary, if these fluctuations are un- 
correlated, coherent EET cannot compete with decoherence 
and the only mechanism left for enhancement of EET seems 
to be dephasing as previously studied in Isj-Q-Sl- The ex- 
perimental results show that for chromophores connected by 
the conjugated polymer backbone there are indeed correlated 
fluctuations in the electronic transition frequencies detected 
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by TTAD. As the authors in 11511 point out, these correlations 
imply that coherent coupling of the electronic degrees of free- 
dom of the chromophores with the vibrational motion, charac- 
teristic of the polymer backbone, must play an important role. 
In fact, they conjecture that the common structural framework 
of the macromolecule is the agent that introduces the correla- 
tions that, in turn, preserves the coherence and permits coher- 
ent EET. 

Of particular relevance to this paper is the recent study of 
Perdomo et al.,|[T3]. They consider an engineered semicon- 
ductor quantum dot array and show that carefully controlling 
the spectral density of the phonon environment can both en- 
hance and suppress energy transfer In this paper we show 
how even a single collective vibrational mode of a nanome- 
chanical resonator can be used to achieve the same effect. 
Our model has the advantage that it can, in principle, be 
tested in the laboratory under well-controlled conditions us- 
ing systems such as nanomechanical resonators coupled to 
quantum dots (v^ or nitrogen-vacancy centers in diamond 
1 18[] or even trap ped ions within optical cavities in the weak- 
coupling regime il9ll . This is essentially the quantum simula- 
tor approach ll20ll . but it may also provide a novel technique 
for controlling quantum transport in future engineered multi- 
component quantum systems. We do not wish to suggest that 
this is in any way a good model for complex naturally occur- 
ring EET systems, such as the FMO complex or conjugated 
polymers, although clearly our motivation derives from such 
chromophoric molecules. 

Our paper is then concerned with configurational changes 
in EET which is an old topic in chemical physics 112 ill . Many 
papers have explicitly included the vibrational degrees of free- 
dom of biomolecules and studied their effect on varied phe- 
nomena. However, most of them are concerned with the cou- 
pling to a broad range of different vibration modes through the 
embedding environment. On the other hand, the solid state 
analogue discussed here makes use of the internal vibronic 
structure in the spirit of the study presented in 12211 . 



GaAS nanomechanical resonator. As shown by Wilson-Rae 
et al. UtIi for the case of a single quantum dot, the vibra- 
tional motion of the beam can indeed change the excitonic 
energy transition frequency. In the case of a coupled array of 
quantum dots, all coupled to the same vibrational mode, the 
modulation in the energy levels will be highly correlated, but 
at the same time fluctuations in the resonator will induce de- 
phasing on the coherent hopping spoiling EET as discussed 
in [22]. We will allow the nanomechanical resonator to be 
classically driven and thus describe it using an underdamped 
simple harmonic oscillator at non zero temperature. Again we 
stress that this is not intended to be a model for natural molec- 
ular EET systems, which are unlikely to be coherently driven; 
rather they are incoherently excited by thermal motion. How- 
ever, our model does necessarily include correlated fluctua- 
tions in the electronic system due to vibrational motion which 
is a central idea for EET emerging from studies of molecular 
systems. 




Figure 1 : One possible realization of the model in which an array of 
excitonic quantum dots are fabricated inside a a GaAs nanomechan- 
ical resonator. The resonator is driven on resonance by a classical 
driving field, and is damped. At very low temperatures, the beam is 
driven to a steady state which is close to a coherent state. The strain 
induced by the deformation of the nanomechanical resonator shifts 
the energies of each dot ITtIi . The corresponding physical parame- 
ters are given in the appendix. 



MODEL 



We start the discussion of our model by considering a typ- 
ical setup used to study the transport of energy between two- 
level dipole systems. It consists of a network of N sites that 
can exchange excitations by dipole-dipole coupUng iIsI-ItIS 



N 



H 



N 



J = l 



(1) 



We now suppose that the dipoles are coupled to a single vibra- 
tional mode. This might be realized for example by coupling 
quantum dots to a single flexural mode of a nanomechanical 
resonator lEsll in such a way that the vibration can modulate 
the on-site energies of the quantum dots. A specific model of 
such a system is shown in Fig.© where a linear array of cou- 
pled excitonic quantum dots are situated in a doubly clamped 



The dynamics of the model, including dissipation, is de- 
scribed by a master equation for the total density operator ( 
(quantum dots plus nanomechanical resonator) l24 1 



dt 



-i[HNV, C] + l{n + l)V[a]Q + -inV[o)]C, (2) 



with 



Hnv = Hn + va)d + e{a^e + i 



N 



where Hj^i given by Eq.([T]i, e is the amplitude of driving 
force on the vibrational mode, and gj is the coupling constant 
between the site j and the motional displacement operator 
q = qo{a^ + a). This coupling generates coiTelated frequency 
fluctuations on the sites as the nanomechanical oscillator vi- 
brates. A very interesting proposal for generating a non driven 
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version of Hmv in the context of trapped ions has recently ap- 
peared 12511 . In Eq.© there is also a term that describes the 
vibrational damping at rate 7, where n is the mean thermal oc- 
cupation of the environment seen by the oscillator, and 2? [a] is 
a super-operator defined as 'D[A]B = ABA^ - {A^ A, B}/2. 
We are using the quantum optics master equation, Eq.Q, 
which assumes that the mechanical motion is underdamped. 
This has been used successfully to describe recent experi- 
ments in nanomechanics ll26ll at low temperature. In the much 
more complicated case of molecular systems there are many 
vibrational modes and a quite different approach is required. 

Much like dephasing enhanced transfer Isf-Q-I^l^ vibration 
enhanced transfer does not occur for all possible Hamiltonians 
of the type given by Q, i.e. for arbitrary LUj, Xjk, and gj. In 
order to present an example of vibration enhanced transfer, 
we consider a simple, yet illustrative, case of a linear chain 
with only nearest neighbor interactions with equal coupling 
constants. In this case, the system Hamiltonian reads 



N 



N 



+e{a^e-"'' + ae""*) + qa{a) + a)Y.,. (4) 

with a^^"^ = and the collective operator = X^jLi 9j^i- 
The vibrational motion of the nanomechanical resonator is 
damped and the driving will settle it into a stable steady state 
oscillation at the driving frequency. For a very low tempera- 
ture of the cooled resonator, this steady state will be close to 
a coherent state, but subject to zero point fluctuations. These 
fluctuations act as a dephasing environment for the quantum 
dots. At finite temperature the steady state is a displaced ther- 
mal state and the dephasing environment is thermal. 

We transform the system Hamiltonian to an interaction pic- 
ture at the oscillator frequency and displace the bosonic op- 
erators a,a^ to subtract the coherent steady state values they 
would have in the absence of coupling to the quantum dots, 
enabling us to treat the steady state of the vibrational motion 
as purely thermal. To this end, we first apply the displace- 
ment operator D{f3) = exp(/3a^ — /3*a), with /3 = — i/3o = 
— 2ie/7, such that the master equation (O will assume the 
form [p = D{-P)ul{t)CUi{t)Dm 



dp 
dt 



C2P + CcP + ClP, 



(5) 



where 



C2P 



Cip 



N 



N 



~2go/3o sini^t 



= -f{n + l)V[a]p + -/nV[a'']p 



(6) 
(7) 
(8) 



In physical terms \P\qo is of the same order as the max- 
imum coherent displacement of the driven nanomechanical 



resonator. Equivalently is the average phonon number 
excitation of the nanomechanical resonator in the steady state. 

In the case of EET in molecular systems, the electronic de- 
grees of freedom are coupled to many vibrational modes. In 
that case, the description usually proceeds by treating the vi- 
brational motion as some kind of heat bath leading to dephas- 
ing of the electronic excitations. In the simplest model this is 
treated as if each site was coupled to an independent heat bath 
II5IQI27, 28], but in recent years it has become clear that this 
is not a good model for real light harvesting systems and tem- 
poral and spatial correlations in the dephasing environments 
arise due to the common vibrational motion of the protein 
cage in which the chromophores are situated IT3 , 15 , U 32 ] . 
In the simpler physical context of this paper there is only one 
vibrational degree of freedom, which can hardly be treated as 
a heat bath. However we can make contact with other mod- 
els by assuming that the vibrational motion is rapidly damped 
on a time scale of the dipole-dipole coupling and the dipole- 
vibrational coupling. 

If the vibrational degree of freedom relaxes rapidly to a 
steady state on the time scale for energy transfer through 
dipole-dipole coupling (A~^), we can adiabatically eliminate 
this motion 133113411 . Notice that if = for all j, the state of 
the oscillator would be a coherently displaced thermal state in 
the limit t >> 7"^ (equilibrium state). In the adiabatic limit, 
the excitonic systems will thus respond as if it were coupled 
to a heath bath; however we need to take care not to average 
away the coherent driving of the mechanical motion in our 
model. We use a version of the Zwanzig projection operator 
method due to Gardiner and Eschmann 13511 . The Zwanzig 
projection operator method then gives the following equation 
for the quantum dots alone. 



&^C2(T+ / dt'TTi[Ccit)exp{Cit')Ccit-t'){pss(E)CF)]. 

Jo 

(9) 

where pss is the steady state of the vibrational mode in the 
absence of the interaction with the quantum dots. This last 
expression allows us to find the master equation for the quan- 
tum dot system in terms of timecorrelation functions of the 
vibrational system evolving under Ci. We find here that 



cr] +4rX'[I]^]cr, 



(10) 



where T = q'^{2n + l)/7 and He is the effective Hamiltonian 
for the network of quantum dots which is given by 



N 



TV 



aEK'^-^'+<^''^-). (11) 



with the time-dependent coefficients Xj{i) ~ '^j/^ — 
2gjf3oqo sin lyt. We can see that the coupling with a driven and 
heavily damped oscillator not only introduces frequency mod- 
ulations but also collective dephasing ll36ll . The role played by 
the external drive on the nanomechanical resonator is now ev- 
ident. It allows the control of Xj{t) through the parameter 
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/?o which is a function of the externally controlled amplitude 
driving e and the damping 7. Our goal in next section is to 
show that by varying /?o, it is possible to find more favorable 
environments for EET when compared to the case of no drive, 
/3o = 0. All simulations in this paper will be performed using 
Eq.dinil. 



RESULTS 

In order to find how efficiently an excitation injected at one 
end of the chain (donor) reaches the other end (acceptor) we 
include an irreversible mechanism of energy loss (sink) in the 
acceptor at the site j = N. By considering only conditional 
evolutions where no absorptions occurred up to time t, we can 
evaluate the rate of energy transfer as the emission rate into the 
sink R{t) = 2K(*c(i)|cr+cr^|«'c(0) S ■ This can be easily 
done by adding a non Hermitian term —iK\e)N{e\ to Hamilto- 
nian (fTTT i and evaluating the time evolved state |^'c(t)) where 
the index c reminds us that this is evolution under the con- 
dition of non-absorption. From this rate, the efficiency for a 
fixed time interval t is then calculated as 



(12) 



Efficiency = / R{t )dt . 



Now we define the basis B = {\ggg--.g) = |0), \egg...g) = 
\i),\geg-g) = \2),\gge...g) ...\ggg...e) = '\N)} 

which is clearly invariant under the action of Hg. The density 



operator in this basis is then written as <t = Sfj=o '^^j I*) 01' 
and the set of coupled differential equations for its matrix ele- 
ments aij obtained from ( fTOl i is given by 

&ij = {-2i{xi{t) - Xj(t)) + 4rGy - k[(5jjv + ^iN^Oij 

— FjNCTij+l], 



in the case i, j ^ where Fij = 1 — 5ij, and 

+i\FjiaQj^i + i\FjN(7aj+i, 



(13) 



(14) 



for the elements associated with the state |0), where Gjj = 
-2{g,- gjf with .90 = 0. 

Using the above equations we numerically obtain the effi- 
ciency for = 6 quantum dots using arbitrary units m = 
h = V = 1 (see the appendix for the physical parameters 
pertinent to a physical example). Please notice that this leads 
to 90 = 1/a/2- In Fig.®, we can clearly see that there are 
amplitudes /3o leading to an improvement on the efficiency 
when compared to the case with no oscillator [constant curve 
(solid)]. In order to analyze the robustness of this mechanism 
of EET, we have also performed a couple of simulations con- 
sidering either Uj or gj as random variables. First, we con- 
sider the same values of gj used in FiglH but now the quan- 
tum dots' frequencies are chosen from independent Gaussian 
distributions. In Fig.©, we show the averaged efficiency over 





1.0 




0.8 








a 


0.6 


2 




'0 

£6 


0.4 


H 






0.2 




0.0 



0.0 0.5 1.0 1.5 2.0 2.5 

Figure 2: Efficiency as a function of /3o for an integration time t — 
3000A. The network frequencies are — 1.0 and cji = a;2 = 
ijJ4 = 1^5 = cijs = 0, the couplings between quantum dots and the 
nanomechanical resonator = 1.5 and 31 = 32 = (74 = Ss = 
ga ~ 0.5, decay constant to the sink k = 0.2, mean number of 
thermal phonons n = 5, and inter-quantum dots coupling A — 0.1. 
The different curves correspond to 7 equals to 1.1 x 10^ (dashed), 
1.1 X 10^ (dotted), and 5.5 x 10^ (dot-dashed). 




Figure 3: Averaged efficiency over 1000 realizations as a function of 
Po for an integration time t — 3000A. The network random frequen- 
cies are chosen from six independent Gaussian distributions having 
mean values ojg — 1.0 and coi — UU2 = ijJ4 = 1^5 — cde ~ 0, 
and standard deviation 0.1, respectively. Decay constant 7 is cho- 
sen as 7 = 1.1 X 10'', and the couplings between quantum dots and 
nanomechanical resonator as well as other parameters are the same 
as in Fig.©. 



1000 realizations. When compared to the Fig.©, one can see 
that the effect of this randomness on the frequencies is to de- 
crease the averaged efficiency. It is interesting though that 
there are still values of /3o 7^ leading to an improvement of 
efficiency. Finally, we now maintain u>j fixed as in Fig.©, but 
we choose the couplings as random variables again according 
to independent Gaussian Distributions. The result is shown 
in Fig.®. The randomness not only decreases the maximum 
efficiency but also tends to restrict the values of /?o leading to 
an enhancement of EET by narrowing the peak. Nonetheless, 
there are still many values of /?o for which EET is favored. 
These results indicate that vibration enhanced EET is effec- 
tive even in the presence of slight randomness. 

If we trace out the first five quantum dots, the resulting den- 
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Figure 4: Averaged efficiency over 1000 realizations as a function 
of Po for an integration time t — 3000A. The couplings between 
the quantum dots and the nanomechanical resonator are chosen from 
six independent Gaussian distributions having mean values gs — 1.5 
and gi ~ g2 = g4 ~ gs = ge = 0.5, and standard deviation 0.3, 
respectively. Decay constant 7 is chosen as 7 = 1.1 x lO'', and the 
quantum dots' frequencies and other parameters are the same as in 
Fig.©. 



sity matrix of the sixth quantum dot will have the off-diagonal 
element given by aoe {t). In order to investigate whether or not 
the coupling to the vibrational degree of freedom can partially 
preserve the coherence of EET, we have calculated | uoe (0 1 for 
the initial state \(f){0)) = (|e)i + \g)i)\ggggg)23456/ V^- The 
result is shown in Fig.©. It can be seen that, for short times, 
the presence of the oscillator (dotted) does indeed improve 
the coherence when compared to the case with no oscillator 
(solid). In this plot, we have kept the sink coupled to the last 
quantum dot k ^ 0. Otherwise, the coherence would assume 
much higher values. 
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Figure 5: Time evolution of coherence for the initial state 10(0)) = 
+ \g)i)\ggggg)23i5e/V2- The network frequencies are oja — 
1.0 and LJi = LJ2 = i.tJ4 = (XJ5 = (xJ6 = 0, couplings between 
quantum dots and vibration mode gi ~ g2 = gi = g4 ~ ga ~ ge ~ 
0.0 (solid), i.e. absence of vibration, gs = 1.5 and gi = g2 = g4 ~ 
Ss = 56 = 0.5 (dotted). The decay constant to the sink is k = 0.2, 
decay constant of the oscillator 7 = 1100, fio = 0.65, mean number 
of thermal phonons in the environment n = 5, and inter-quantum 
dots coupling A = 0.1. 

This is an example of enhancement of quantum transport 
in a network of two-level systems by coherent coupling with 



a driven resonator. The role of incoherent vibronic coupling 
(phonon-induced decoherence) in EET has already been in- 
vestigated H, 22]. Our simplified model, valid in the inter- 
mediate coupling, has shown the possibility of coherent EET 
with partial preservation of phase information, something un- 
likely to take place in noise assisted EET. 



DISCUSSION AND CONCLUSIONS 

As one can see from Fig|2] energy transport occurs at par- 
ticular values of the mechanical resonator driving strength. To 
understand this we neglect the decoherence due to the ground 
state fluctuations of the vibrational mode as well as absorption 
at the special site and consider only the network Hamiltonian 
which we write as 



N 



N 



where the modulated site energies are defined as 

UJj{t) = ujj — AgjqoPQ sinz^i. 
We now transform to an interaction picture via 



Uoit) = exp 



■ N t 
j=i "'0 



(15) 



(16) 



(17) 



to obtain the interaction Hamiltonian 



N 



(18) 



where 9{t) = Aujjt — 4iAgj{f3oqo/v) {cosft — 1), Aujj = 
Ljj — and Agj = gj — gj+i- Expanding the oscillating 
exponential as a Fourier series using the generating function 
for Bessel functions J„(x), we obtain 

n— 00 N 

Hj{t) = X J2 ^[J'«e'(^"^-"''>V>i+V/i.c.], (19) 

n— — 00 j—1 

where we defined the coefficients 

Jn - (i)-"J„(4A5,/3ogo/j^)e^^''^^^"««/". (20) 

It is now apparent that there will be resonances when 

Aujj ~ nv. (21) 

We thus expect that when the differences between site fre- 
quencies is of the same order as the mechanical vibrational 
frequency, the resonant terms in Eq.(fT9b will dominate the in- 
teraction, provided the Bessel function coefficients have suffi- 
cient weight at those frequencies. These resonances are essen- 
tially due to the sidebands written on the electronic transitions 
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by the phase modulation implicit in Eq.lfTSIl. The Bessel func- 
tions are a function of the driving amplitudes, /3o, and also 
the coupling inhomogeneities, Agj. Note that if the system 
is completely homogeneous so that Aujj = and Agj = 0, 
there are no resonances. Inhomogeneity plays a crucial role. 
Given some inhomogeneity, the resonances can be controlled 
by the external driving strength through the dependence on 
/3o. 

We can see the effect of the Bessel function weighting on 
the resonances in Fig |2] In the dimensionless units of this 
example the vibrational frequency is = 1 while the largest 
difference in site frequencies is Aaj2 = — Awa ~ 1.0 so a 
first order resonance is possible. The first two zeros of the 
first order Bessel function occur at xi — 3.83 and X2 — 7.02. 
This would imply a suppression of the resonance for the corre- 
sponding values of /?o ~ 135 and /3o ~ 2.5 and we do indeed 
see a suppression of EET near these values. Note that it is 
the product of the driving strength /3o and the coupling differ- 
ences Ag^ that matters. If the coupling is weak, one needs 
higher driving strengths to see the same phenomenon. This is 
illustrated in Fig.®, where the zeros of the first order Bessel 
function would now imply in suppression of EET at /3o ~ 45 
and /3o ~ 83, in very good agreement with the numerical plot- 
ted curve. Likewise, larger vibrational frequencies require a 
higher driving strength. 
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Figure 6: Same as Fig.([2l( but now with 33 = 0.03 and gi ~ g2 ~ 
gi = gs ~ ge = and equal to 7 — 1.1 x 10^. This value of g is 
consistent with fivll . as discussed in the appendix. 



To summarize, we have shown that for an inhomogeneous 
linear chain of coupled quantum dots in a harmonically driven 
vibrational structure, particular driving strengths enhance the 
quantum transport of excitations across the network. Future 
work may cover the inclusion of additional vibrational degrees 
of freedom and move beyond nearest neighbor coupling. We 
have also shown that vibration enhanced quantum transport 
can partially preserve quantum coherences during electronic 
energy transfer. This fact may motivate future research on 
propagation of quantum information in this system. 
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APPENDIX 

In this appendix we calculate the model parameters for a 
specific realization based on quantum dots in GaAs, following 
the proposal of Wilson-Rae et al. IitIi . The coupling between 
the vibrational mode and a single quantum dot in their model 
is given by 



VD 



hvr] 



(22) 



where 77 is a dimensionless parameter given by Eq. (3) in UtIi . 
which we also take as 77 = 0.06. Comparing this to Eq.© in 
our model we see that 



n-9 = TT- 



(23) 



where qo is the rms position in an oscillator ground state and 
is given by 



go = 



2mi' 



1/2 



(24) 



s ^ and 50 ^ 5 X 10 ^'^ m. This 
^ in standard units. In the units 



We use the results from IitIi for a GaAs beam of length, L ~ 
1 jim, width 6 = 85 nm and depth d ~ 30 nm, the mass is 
of the order of m ^ 1.4 x 10~^^ kg, the fundamental mode 
frequency is ;^ = 1.2 x 10 
implies that g 36 x 10® s 
of this paper such value of g translates as g = 'j]/2 = 0.03, 
and this was considered in Fig.©. 

The damping rate for the nanomechanical resonator is de- 
termined by the quality factor, Q, of this resonance. If we 
take the rather conservative value of Q = 10^, then at a res- 
onance frequency quoted by Wilson-Rae et al. lflTll . we have 
that 7 ~ 10^ s~^ leading to gqo/j ^ 1 which is the weak 
coupling regime. In this regime, the oscillator is heavily 
damped compared the coherent coupling with the quantum 
dots. Should this limit not be sufficient to safely eliminate 
the vibrational motion adiabatically, a lower mechanical reso- 
nance frequency or a lower Q resonator may be used. 

The site energy for a GaAs quantum dot is of the order of 
hiOi ^ 1 meV 1 38]. From Perdomo et al. flill we take the ratio 
of the coupling strength to the excitation energy of the quan- 
tum dot to be to be of the order of — ~ 0.1 . This corresponds 
to quantum dots separated by about 10 nm. In order to observe 
the kind of resonance discussed in this paper we need the site 
energy inhomogeneity, A (/to;) to be of the order of the vibra- 
tional energy quanta hiy. For the case of = 1.2 x lO^s^^, 
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this would imply Au ^ 1 /le V. This is somewhat smaller than 
typical inhomogeneities that occur due to growth conditions. 
However external electric fields can be used to tune quantum 
dots to a considerable degree js^, so it may be possible to 
set up conditions to observe the phenomenon discussed in this 
paper using coupled InAS/GaAs quantum dots in a nanome- 
chanical resonator. A detailed description and theory of such 
an experiment will be given in a future paper. 

Note added in proof: While being considered for publica- 
tion, some independent studies related to this work appeared 
in the Uterature |40] 
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